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We examine a data-rale synchronization system capable of operating 
in two modes: (i) in master-to-slave mode when the data stations are con- 
nected by digital transmission facilities, and (ii) in slave-to-slave mode 
when the data stations are connected by analog transmission facilities. 
The first part of this paper determines the steady-state behavior and the 
transient response in the master-to-slave mode. The results show that the 
system is well behaved in the transient stage, and that the steady-state 
behavior is satisfactory. From the transient analysis, the buffer size require- 
ments of the system and the counter size requirements of the rate-locked 
loops are determined. Formulas are developed from which the start-up 
time of the system can be estimated. 

The second part of this paper examines the behavior of the system in 
the slave-to-slave mode. It is shown that the data stations can settle to the 
same steady-stale signaling rate, and this signaling rate is determined. 
The dependence of this signaling rate on other system parameters is ex- 
amined. It is shown that the system can be easily designed such that the 
steady-state signaling rate will lie within desired limits. (This is so regard- 
less of the starting sequence, the initial system conditions, and time delays 
in the communication channels.) 

I. INTRODUCTION 

When data stations arc connected by wholly digital transmission 
facilities, it is most efficient to slave the clocks at the data stations 
to a master clock. To perform this operation, hereafter referred to as 
master-to-slave operation, an interface unit at the data station extracts 
timing pulses from the incoming data stream. These timing pulses 
are passed through a phase-locked loop to eliminate noise and jitter. 
The output of the phase-locked loop controls the signaling rate of the 
data station. 
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Unfortunately, a technical problem arises when data stations are 
synchronized in the above manner. Before digital systems evolve into 
a well-connected network, data stations are also often connected by 
wholly analog transmission facilities. When two data stations equipped 
to operate in the master-to-slave mode are connected by analog facilities, 
each station will regard the clock at the other station as the master 
clock, and the two stations will attempt to mutually synchronize each 
other. This mode of synchronization can be called "slave-to-slave." 
Conventional phase-locked loops 1 which perform well in the master- 
to-slave mode may not perform well in the slave-to-slave situation, 
being unusually sensitive to path-length delays and other system 
parameters. This technical problem can be solved by avoiding the 
slave-to-slave situation in the following manner: 

(i) Informing the data stations when analog transmission facilities 
are used. This will permit the stations to break up the slaving 
paths in the data sets and use their own clocks as the timing 
source. 
(it) Providing a looped connection within the analog system con- 
taining appropriate buffers, and a clock of sufficient accuracy 
to serve as the master for the data stations. 

Unfortunately, these schemes reduce the economic attractiveness 
of the system. Consequently, there is a need for a synchronization 
scheme capable of operating in both the master-to-slave and the 
slave-to-slave modes. 

We analyze a synchronization system which employs digital rate- 
locked loops to determine if it can operate successfully in both modes. 
The phase detector in the rate-locked loop is a multistage counter that 
counts the difference between the number of zero crossings of the input 
signals. Because of this nonlinear counting process, the operation of 
the synchronization system is determined by nonlinear differential- 
integral equations. Such equations do not appear in earlier synchroniza- 
tion studies 2-4 which considered different phase detectors. As will be 
shown, a digital rate-locked loop locks to neither the phase nor the 
frequency of the timing signal, but to the zero-crossing intervals. This 
difference complicates the analysis. We have examined the problem 
without making linear approximations. 5-8 In a previous paper, 9 we 
analyzed, in a rigorous fashion, the steady-state behavior of the system 
in the master-to-slave mode and proved that, in the absence of filtering 
in the rate-locked loop, the slave oscillator will lock to the master 
oscillator exactly. In this paper, it is proved rigorously that if the 
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filter in the loop satisfies a simple condition, the system will reach 
equilibrium (that is, data stations cannot add or delete bits from a 
customer's data stream). Based on this, it is demonstrated that in 
the presence of RC filtering in the loop, the slave oscillator will lock 
to the master oscillator exactly. Following these analyses, this paper 
determines the transient response of the system in the master-to-slave 
mode, and examines the behavior of the system in the slave-to-slave 
mode. Sections II and III of this paper examine the master-to-slave 
mode. Transient response, buffer-size requirements of the system, and 
counter-size requirements of the rate-locked loop are determined. 
Section IV considers the slave-to-slave mode. Steady-state signaling 
rate of the system is determined, with its dependence on the other 
system parameters examined. A simple method of designing the system 
to ensure satisfactory steady-state signaling rate is presented. Section V 
summarizes the results of this paper and may be read next. 

II. MATHEMATICAL MODEL 

In this and the following two sections, we examine the master-to- 
slave mode. Consider two communication stations as depicted in 
Fig. 1. Station 1 (with slave clock) represents a data station. Station 2 
(with master clock) represents a station in the digital transmission 
facility. The master clock at Station 2 emits a timing signal which 
controls the transmission of data from Station 2 to Station 1 (for 
example, Station 2 transmits a digit to Station 1 at every second zero 
crossing of this timing signal). Station 2 transmits to Station 1 at some 
standard rate, say, / 2 digits per second. 

Station 1 receives data from Station 2, and derives from the received 
data a timing signal s 2 (t) = sin (w 2 t + 8 2 ), where o> 2 = 2ir/ 2 and 2 is 
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Fig. 1 — Master-to-slave operation, block diagram. 
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an arbitrary phase angle. The signal s 2 (t) and the output s, (*) of a local 
oscillator are compared in a digital phase detector (Fig. 2). The digital- 
phase detector is a counter which counts the zero crossings of s 2 (t) 
and s,(0, and produces an output proportional to the difference between 
these two counts. Mathematically, this operation can be specified as 
follows. Let it be assumed that the digital phase detector is activated 
at t = 0. Let Ni(t) and N 2 (t) be, respectively, the number of zero 
crossings (both upward and downward zero crossings) of s v (t) and 
s 2 (t) in the time interval to t; then the output of the digital phase 
detector is 



«i(0 = e 1 [N a (t) - #,(0] 



(1) 



where ei is a positive constant (volts/count) and may be called the 
gain of the counter. As depicted in Fig. 2, Ui(t) is passed through a 
filter, and the filter output v t (t) controls the frequency of a voltage- 
controlled oscillator (VCOi). Let o>i = 2irf, be the free-running radian 
frequency of VCOi , then the output of VCOi is 



8,(0 = sin &»!< + «i J »i(r) dr + M 



(2) 



where a x is the gain of VCOi (radians/volt X second). The signal 8 t (t) 
is used to control the transmission of data from Station 1 to Station 2 
(for example, Station 1 transmits a digit to Station 2 at every 2nd 
zero crossing of s,(0). Note that t in (2) represents the phase of s^t) 



S,(t) 



s 2 (t) 




Fig. 2 — Digital phase detector and the rate-locked loop at Station 1. 
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Fig. 3— Illustration of S 2 (l), N t (t), and fa(t). 

at t = 0. Without loss of generality, we may assume ^ d x ^ ir, and 
S e 2 g, t. 

Let us derive an analytic expression for the number of zero crossings 
of s 2 (t) from t = to a particular time instant t. As illustrated in Fig. 3, 
t' is the time instant at which the last zero crossing prior to t = takes 
place, and t" is the time instant at which the last zero crossing prior 
to t = t takes place. It is obvious from Fig. 3 that the number of zero 
crossings in the time interval to t is 



N 2 (t) = 



co 2 r + e 2 



(3) 



Note from (3) that the phase cumulated from t" to t does not con- 
tribute to the value of N 2 (t). This residual phase (or phase quantiza- 
tion error) will be designated ^ 2 (t), i-e., 



^ 2 (t) = w 2 t — i0 2 t' 



(4) 



Equations (3) and (4) hold for all t > 0; therefore, we can replace their 
t by the time variable t. The variation of ip 2 (t) with t is illustrated in 
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Fig. 3. Note that \J/ 2 (t) increases from to ir. When \f/ 2 (t) reaches tt 
radians, a zero crossing takes place; and fait) drops to zero and increases 
from zero again. Clearly, ^ $ 2 {t) ^ *"• Since s 2 (t) is a pure sine wave, 
\f/ 2 (t) is a sawtooth wave. 
From (3) and (4), we have 



N 2 {t) = 



C0 2 * + 02 - M) 



Similarly, one can write the number of zero crossings of 8 t (t) as 

; 1 i + a 1 f vMdr+ d, - <M«) 
Jo 



N^t) = 



(5) 



(6) 



where ^(i) is the residual phase as illustrated in Fig. 4. As can be 
seen, ^ ypi(i) ^ ir. Note that \f/i(t) is not shown as a sawtooth wave 
in Fig. 4 because Si(t) is not a pure sine wave in the transient stage 
after t = 0. 

In this paper, the filter F(s) in Fig. 2 is assumed to be the usual RC 
filter (Fig. 5). Thus, its transfer function F(s) is 1/(1 + sCR). Sub- 
stituting (5) and (6) into (1), and rearranging the equation, we obtain 



Si(t)=SIN[a)jt + ai /vi(T)d7-t0,] 




Z^, 



Fig. 4— Illustration of Si(t), Ni(t), and fr(«). 
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Fig. 5 — RC filter in the rate-locked loop. 



t (t) = kht - a, f' V,{r) dr + 2 - 0, + «M0 - *i(0J 



where 



»-* 



5 = w 2 — Wi . 



(7) 



(8) 



(9) 



We shall use one-sided Laplace transform in the analysis (the words 
one-sided will be omitted). As usual, the Laplace transform of a time 
function will be consistently denoted by the appropriate capital letter. 
For instance, Ui(s) will denote the Laplace transform of u t (t). The 
symbol L[f(t)] denotes the Laplace transform of f(t), and the symbol 
L _I [F(s)] denotes the inverse Laplace transform of F(s). Taking the 
Laplace transform of (7), we obtain 



[_S S 5 o _J 



(10) 



Multiplying both sides of (10) by F(s), using F(s)C/,(s) = 7i(s), and 
rewriting the resulting equation in time domain we obtain 



Vl (t) = L"{h( S ) s 4] + L-{#00f] - L-{#( S )f] 
- L-[ff («)¥,(«)] + L"W.(^1. * > 



(ID 



where 



H(s) = 



ks 



CR(s + r^s + r 2 ) 
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The two roots r t and r 2 are real numbers when 1 — 4CRka t > 0, and 
are complex numbers when 1 — ^CRka t < 0. It can be shown that in 
the second case the frequency of s t (t) overshoots that of s 2 (t) before 
it finally settles. Such an overshooting should be avoided because 
Station 1 may be required to operate in the slave-to-slave mode (see 
Section V). Therefore, throughout this paper we assume 

1 - iCRka, > 0. (12) 

III. STEADY-STATE AND TRANSIENT ANALYSES 

In the master-to-slave mode, we have to consider the following 
questions: 

(i) Can the signaling rate of Station 1 lock to that of Station 2 
in the presence of phase quantization errors? 

(ii) What is the steady-state frequency of VCOi ? 

(Hi) During the transient stage after t = 0, the signaling rate of 
Station 1 can be higher than that of Station 2. Therefore, data 
can be transmitted from Station 1 to Station 2 faster than it 
can be transmitted out of Station 2. Consequently, a buffer 
storage is required at Station 2. What should be the size of this 
buffer? 

(iv) The digital phase detector is a counter that counts the dif- 
ference between the number of zero crossings of s t (t) and s 2 (t). 
How many stages are required in the counter to avoid overflow 
(i.e., to ensure pulling in)? 

We shall first determine the transient response of the system in 
Section 3.1, and then consider these questions in Sections 3.2 to 3.4. 

3.1 Transient Response and Settling Time 
We evaluate the first three inverse transforms in (11) to obtain 

L- l [ff0pJ = CRrMrl - r 2 ) 



[fri - r 2 ) + r 2 e-" 1 - r^' 1 ] (13) 
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L-{H(.)f]- l -[H«t] = fc^^'--n M 

The fourth inverse transform L _1 [//(s)¥ 2 (s)] is more difficult to evaluate. 
We obtain after lengthy manipulations 

ITlffd)*^')] = a(t) + p(t) (15) 

where 

+ , k \-* +e 2 - -^r r V r ". 

Vl - 4Cfl/ca, L r 2 l _ e r ' T J 



(16) 



p(0 = a periodic function of period T, identical with p (t) in the time 
period ^ t ^ T (17) 

+ . M^- — [1 - (1 - e" T M< - OK"'- (18) 

Vl - 4C#/ca, (1 - e r ' T ) 

Note that a(t) is the sum of two decaying exponential terms. When t 
increases, a(t) approaches zero, and only the periodic steady-state 
response p(t) remains. It can be shown that p(t) has zero mean. 

Now consider the last inverse transform L _1 [#(s)^,(s)] in (11). 
Since *,(s) is the Laplace transform of M), L~'[ff («)¥,(«)] can be 
evaluated if ^,(0 can be determined. As illustrated in Fig. 4, ft(Q 
depends on the positions of the zero crossings of s,(/). Furthermore, 
because we are dealing with a closed-loop control system, ^i(<) and the 
phase of s, (t) must satisfy the integral equation (7). In order to determine 
^,(0, one must simultaneously consider (7) and the zero crossings of 
8i(t). The mathematical problem is extremely complex and it is im- 
possible to obtain a closed-form expression of ^(0 for all t. Con- 
sequently, the inverse transform L _, [#(s)*i(s)] cannot be evaluated 
in closed form. However, we have obtained a tight upper and lower 
bound for its value as follows: 

|L- 1 [H(s)^ 1 (s)]| < 2 6l . (19) 
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We have obtained the closed-form expression of the first four com- 
ponents of Vi(t), and tightly bounded the fifth component of Vi(t). This 
gives the transient response of the system. Note from (13) to (19) 
that transients in v } (t) either decay exponentially or can be bounded 
by a small number. Thus, the system is well behaved in the transient 
stage. Furthermore, from these equations, one can plot Vi(t) vs t, and 
easily estimate the settling time of VCOi . The settling time of VCOi 
can be rather long when CR is large. For example, consider the first 
term L~ 1 [H(s)5/s 2 ] in v x {t) (this is usually the dominating term in 
Vi(0). From (13) it can be rewritten as 

Since r t > r 2 > 0, the last term in the right-side bracket is negative 
for all t. Thus, the convergence of l~ 1 [H(s)8/s 2 ] is even slower than the 
convergence of the time function 1 — e~ r ''. This clearly shows that 
Vi(t) converges slowly when the filter time constant CR is large. 

3.2 Steady-State Frequency of VC0 1 

Now we answer the first two questions at the beginning of Section III. 
First, we have found that the signaling rate of Station 1 will lock to 
that of Station 2 in the presence of phase quantization errors. The 
proof of this is complicated, and is given in the Appendix. In this 
section, we examine the steady-state frequency of VCOi , and point 
out an important difference between digital and analog phase detectors. 
The instantaneous frequency of VCO, is [ui + ot l v 1 {t)/2ir\. In order 
to see if it approaches a fixed steady-state value, we evaluate 
lim,-,,, «i(0, which can be found by evaluating the limits of the five 
inverse transforms in (11). As shown in the Appendix, when signaling 
rate of Station 1 locks to that of Station 2, the zero crossings of Si(t) 
and s 2 (t) will alternate with probability one, and yf/^t) will be a periodic 
function of period T. This means that L~ 1 [H(s)\f/ 1 (s)] also approaches 
a periodic function of period T. Let this periodic function be denoted 
by q(t). Then, one can show from (11) that 

/Instantaneous fre-\ a t ... . a, ... , nn . 

Iquency of VCO, ) = U ~ 2* *® + Y« q{t) ' (20) 

When signaling rate of Station 1 locks to that of Station 2, zero 
crossings of s,(0 and s 2 (t) alternate. Therefore, fa(t) 5* yj/ 2 (t) and 
P(t) ** ?(0- Thus, from (20), instantaneous frequency of VCO x does 
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not lock to the master clock frequency f 2 . Instead, it is a periodic 
function f 2 - (aj2ir)p{t) + (ctj2ir)q{t). The output s l (t) of VCO, is 
a periodic wave with the same period as s 2 (t); however, it is not a pure 
sine wave as one would expect from experience with analog-phase 
lock loops. The digital loop locks to neither the instantaneous frequency 
nor the -phase of the incoming signal s 2 {t). It locks only to the rate of zero 
crossings of s 2 (t). For this reason, it should be referred to as a digital 
rate-locked loop, rather than a digital frequency-locked loop or a digital 
phase-locked loop. This difference between digital and analog loops 
should be noted in the applications. 

3.3 Size of the Data Buffer at Station 2 

As described in Section II, Station 2 transmits to Station 1 at a 
standard rate of f 2 digits per second. In general, Station 1 is also re- 
quired to transmit to Station 2 at this standard rate. To achieve this, 
Station 1 transmits a digit to Station 2 at every second zero crossing 
of Si{t) [this enables station 1 to transmit also at the standard rate 
when Si(t) is synchronized to s 2 (t)]. 

Usually, Station 2 relays the data it receives from Station 1 to another 
station at the standard rate of f 2 digits per second. Thus, when the 
system is in synchronization, data is transmitted to Station 2 at the 
same rate as it is transmitted out of Station 2. However, when Station 1 
is first synchronized (that is, during the transient stage of synchroniza- 
tion), the transmission rate of Station 1 can be higher than f 2 . Con- 
sequently, during the transient stage, data can be transmitted from 
Station 1 to Station 2 faster than it can be transmitted out of Station 2. 
This means a data buffer is required at Station 2. In this section, we 
determine the size of this buffer. 

As defined in Section II, JV x (0 is the number of zero crossings of 
Sx(t) in the time interval to t. Since Station 1 transmits a digit to 
station 2 at every second zero crossing of s,(0, the number of digits 
transmitted from Station 1 to Station 2 in the time interval to t is 
Nx(t)/2 or Ni(t) — 1/2, depending on whether Ni(t) is even or odd. 
To simplify our discussions, we shall use the following definition through- 
out this paper. 

Definition: For any positive number a, (a) denotes the integer im- 
mediately less than a when a is not an integer, (a) = a when a is an 
integer. 

Using this definition, the number of digits transmitted from Station 1 
to Station 2 in the time interval to t is {N x (t)/2). The number of 
digits Station 1 should transmit in this time interval is (N 2 (t)/2). If 
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(Ni(t)/2) is larger than (N 2 (t)/2), a buffer would be required at Station 2 
and the buffer size is (N x (t)/2) — (N 2 (t)/2) digits. It can be shown 
from the previous equations that 



-K + 3] 



< JVt(0 - N 2 (t) < — + 3. (21) 

e x a x 



Since the buffer size is (N x (t)/2) — (N 2 (t)/2), we obtain from (21) 
Buffer Size < ^ \- 2. (22) 

2e x a x 

Equation (22) gives an upper bound for the buffer size. It can also be 
shown that in order to prevent overflow the buffer size must be greater 
than (io'/2e x a x ) — \. Combining this with (22), we have 

- \ < Buffer Size < ^— + 2. (23) 



2e x a x 2e x a x 

Let us define B = (&>' '/2e x a x + 2). It can be seen from (23) that the 
buffer size is J5, B — 1, or B — 2. Thus, the buffer size is determined 
to within two digits. Since the two-digit difference is negligible, one 
may use the simple formula 



Buffer Size = B = l-^- + 2 ) digits. 
\2eitti / 



(24) 



As explained at the end of Section II, in this paper we use the con- 
straint 1 — 4:CRka x > 0. From this constraint, we can rewrite (24) as 

Buffer Size = Z 2 ^^ + 2) digits (25) 

where ft = kCRe^/v. Clearly, < ft < 1. Equation (25) will be 
used in later discussions. 

3.4 Counter Size of the Digital Phase Detector 

The counter in the digital phase detector counts the difference between 
JV 2 (0 and NxQ). Now we determine the counter size so that the counter 
will not overflow when the maximum positive count or negative count 
is reached. Consider first the case of negative counts. It can be shown 
that N x (t) — N 2 (t) can be larger than <J /e x a x . It has been shown in 
the preceding subsection that N x (t) — N 2 (t) must be less than <J le x a x 
+ 3. Thus, if we define 
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the integer N^t) - N 2 (t) can be as large as N - 2, but will not exceed 
N. Therefore, the counter will not overflow if it can count Ni(t) — N 2 (t) 
up to N t (t) - N 2 (t) = N. 

Next, consider positive counts. One can show that the counter will 
not overflow if it can count N 2 (t) - N l (t) up to N 2 (t) - N t (t) = N. 
Combining the two cases, we see that the counter will not overflow if 

Counter Size = ±N counts = ±(-Jj h 3 ) counts (26) 

where fr is defined after (25). 



IV. SLAVE-TO-SLAVE SYNCHRONIZATION USING DIGITAL RATE-LOCKED 
LOOPS 

In this section, we consider mutual synchronization between two 
data stations where each station regards the clock at the other station 
as the master clock. Such a mutual synchronization is usually called 
slave-to-slave synchronization. 

A mathematical model of slave-to-slave synchronization is depicted 
in Fig. 6. The local oscillator at Station 1 (VCO, in Fig. 6) emits a 
timing signal S n (t) which controls the transmission of data from 
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Fig. 6 — Slave-to-slave synchronization with digital rate-locked loops at both stations. 
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Station 1 to Station 2. [For example, Station 1 may transmit a digit 
to Station 2 at every second zero crossing of 8u(t).] Station 2 derives 
from the received data a timing signal S 12 (t) and compares *Si 2 (0 with 
its local oscillator output S 22 (t) at the digital phase detector. The 
digital phase detector is essentially a counter which counts the zero 
crossings of S 12 (t) and S 22 (t) and produces an error signal u 2 (t) pro- 
portional to the difference between these two counts. The error signal 
u 2 (t) is passed through a filter F 2 (s) to control the frequency of VC0 2 . 
Thus, in this fashion, Station 2 adjusts its clock rate toward that of 
Station 1. Similarly, as depicted in Figure 6, Station 1 regards the 
clock at Station 2 as the master clock and adjusts its clock rate toward 
that of Station 2. 

Practically, it is impossible to activate the two counters at the two 
different stations at the same time instant. Therefore, in this study, 
we consider an arbitrary starting sequence as follows: 

(i) At an arbitrary time instant t t , either the counter at Station 1 

or the counter at Station 2 is activated. 
(ii) The other counter is activated at an arbitrary later time instant 
t 2 (t 2 > tt). 

For analytical purpose, we shift the time origin such that t 2 = 0. 
We shall analyze the behavior of the system for t > 0. 

Let o>i be the free-running radian frequency of VCOi , then we can 
write 



s n (t) = sin Pll (t) 

= sin Wit + oil I v^t) dr + #n 



(27) 



and 



s i2 (t) = sin [p n (t - r 12 )] (28) 



where t 12 is the time delay introduced by the channel. Similarly, the 
free-running frequency of VC0 2 is denoted a> 2 and 

S22O = sin p 22 {t) 

= sin co 2 t + a 2 J v 2 {t) dr + B 2i (29) 

and 

Sai(0 = sin [p 22 (t — T 2 i)]. (30) 
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We define N it (f) as the number of counts from s,,(0 in the time 
interval to t. Let us first derive an analytical expression for N n (t). 
From (27), we can write 



«i* + «i [ «i(t) dr + B u = MirN n (t) + iMO 
«/o 



(31) 



where the parameter M is defined by: M equals one when the counters 
at the two stations count both the upward and downward zero crossings; 
M equals two when the counters count only the upward (or downward) 
zero crossings. The last term $ n (t) in the above equation represents 
phase quantization errors and ^ ^ u (0 < Mr. From the above 
equation we have 

NUt) - jjfj [«i* + «. f »M dT + *" - *»(0J ■ (32) 

Similarly, one can write analytical expressions for N l2 (t), N 2l (t), and 
N 22 (t). 

At Station 1, the digital counter output is 

Wi(0 = u,(0) + e,[iV 2l (0 - AT u (0] (33) 

where Wi(0) is the initial count at t = 0. Let us define 

*i = it (34) 

Mir 
6 = oj 2 — o>i (35) 

0, = 21 - 9„ (36) 

fla.(s) = e- r " [ v 2 (t)e-' dt. (37) 

The filters F,(s) and F 2 (s) in Fig. 6 are assumed to be the usual 
J?C filters, i.e., 

'•» - r+fes; (38) 

«* = r+kfi;- (39) 

At Station 2, the digital filter output is 

u 2 (t) = u 2 (0) + e 2 [N»(t) - N 22 (t)] (40) 
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where u 2 (0) is the initial count at t = 0. Furthermore, we define 

fc = ifc <«> 

$2 = $12 — $22 (42) 

BM = e" r " j Vtifie- 1 dt. (43) 

From the previous equations, we can determine Vi{s) and V 2 (s). 
The results are: 



v m _ A ^ . Ms) , A 3 (s) A 4 (s) A 5 (s) . 

Vl{s) ~ B(s) + B(s) + B(s) + B(s) + B(s) (44) 



where 



Ms) = \uMEM + M£M + , i(fflCiBi , i(f) + >*R.WM l 

|_ s s s J 

\_ S S S J 

•[^^ /?,(«)] 
A 4 (s) = [-k 2 * 12 (s)F 2 (s) + fc 2 ^22(s)F 2 (s)][fc 1 a 2 ? -^i5' 1 (s)J 
A 8 (s) - k t 8 £j& [l + fc 2 « 2 ^ (1 - a"'")] 



B(«) = 1 + Jfc iai ^ + A;^ ^ 
s s 

+ *««*«,£&£& [1 -«-<'"«">]. 

S S 

Similarly, we obtain 



v m _ ^«00 . Ms) , A B (s) A«(s) AM ,ar\ 

KiW B(s) + B(s) + B(s) + fi(«) + B(s) (45) 
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where 

. TuMEM + hhEM + , ( „ )CM( ,) + h^&Mm] 

\ S S *> J 

A 7 (s) = [-k 2 *MF 2 (s) + /c 2 ^ 22 (s)F 2 (s)][^l + k^ ^J 

Ms) = \uMlM + M£M + ^W^) + fc ^^> l 

L s s s J 

.[w«M] 

A 10 («) = -fc 2 5 ^ [l + ft*, ^ (1 - e- T -)J 

Note that our problem is not solved. Equation (44) is not a closed- 
form solution of Vi(s) because ¥<*(*), which appears in A 2 (s) and 
A 4 (s), depends on 7,(s) and V 2 (s) (the phase quantization errors ^,,(0 
depends on v t (t) and y 2 (/)). Similarly, (45) is not a closed-form solution 
of V 2 (s). These equations will, however, enable us to examine the 
steady-state behavior of the system in the following sections. 

4.1 A Steady-State Solution o] Signaling Rates 

As described previously, the zero crossings of s u (t) are used to control 
the transmission from Station 1 to Station 2 (for example, Station 1 
may transmit a digit to Station 2 at every second zero crossing of 
s u (t)). Similarly, the zero crossings of s 22 (0 are used to control the 
transmission from Station 2 to Station 1. Therefore, to determine the 
steady-state signaling rates of these two stations, it suffices to determine 
the steady-state distribution of the zero crossings of s n (t) and s 22 (t). 
To facilitate our discussion, let us first introduce the following definition. 
Definition: s (l) denotes a sine wave sin u Q t with 

wo = i r~ ? TT» — T. — 7 Z ^ [wifca«a + t*Jhf*i 

K 1 a i ■+- K 2 a 2 ■+■ /CiO!i«2<*2t r 12 T T 2 \) 
+ aj 1 fc 1 a 1 fc 2 a 2 Ti 2 + w 2 k 1 a 1 k 2 a 2 T 2 i 
+ [Mi(0) + Mi + fcia 2 fl2i(0)]/c 2 a,a 2 
+ [1*2(0) + k 2 d 2 + JbtaiB M (0)]JbiffiaJ. (46) 
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Based on (44) and (45), a steady-state solution of the zero-crossing 
distribution has been obtained. 10 In order to conserve space, let us 
omit the lengthy derivations and write only the results as follows: 
A Steady-State Solution: When the counters at the two stations count 
both the upward and downward zero crossings of s ,-,(£), i t J = 1, 2, 
the upward and downward zero crossings of s,,(i), *, ; = 1, 2, are 
uniformly spaced when t — * <x> and the time interval between each two 
consecutive zero crossings of Sa(t), h i ™ 1, 2, is identical with the 
time interval between each two consecutive zero crossings of s a (t). 

If the counters count only the upward (or downward) zero crossings, 
the above solution should be modified: When the counters at the two 
stations count only the upward (or downward) zero crossings of s,,(Z), 
i, j = 1, 2, the upward (or downward) zero crossings of s,,(£), i, j = 
1, 2, are uniformly spaced when t — > °o , and the time interval between 
each two consecutive upward (or downward) zero crossings of s,,(i), 
i, j = 1, 2, is identical with the time interval between each two con- 
secutive upward (or downward) zero crossings of s (t). 

4.2 Analysis of the Steady-State Signaling Rate 

In this section, we show that the system can be easily designed such 
that the steady-state signaling rate lies within desired limits. 

Before the two stations are mutually synchronized, s n (t) is sin u^t 
and the signaling rate of Station 1 is hui digits/second, (h is a pro- 
portionality constant. For example, h = l/2ir when Station 1 transmits 
a digit at every second zero crossing of s n (t).) Similarly, before the 
two stations are synchronized. s 22 (t) is sin u 2 t and the signaling rate of 
Station 2 is hu> 2 digits/second. When the two stations are mutually 
synchronized, s u (t) and s 22 (t) have the same zero-crossing distribution 
as s (t) = sin co t and the signaling rates of the two stations are ho)„ 
digits/second. The synchronization is satisfactory if hw is sufficiently 
close to hcoi or hu 2 . More specifically, the steady-state signaling rate 
is satisfactory if 

/iwi — e < haio < hui 2 + e (47) 

when to, < w 2 , and if 

hw 2 — e < ho) < hcoi + c (48) 

when o) 2 < wi • The number e is a prescribed small number. 

As can be seen from (46), u> depends on a>i , w 2 , and the following 
parameters: gains e x and e 2 of the two counters, gains a x and a 2 of the 
two oscillators, initial counter outputs w,(0) and u 2 (0), initial phases 
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d l and 2 , initial filter outputs v^t) and v 2 {t), and the time delays t 12 
and r 2 i in the communication channels. Since w depends on so many 
parameters, it is not immediately clear whether o> satisfies the specifica- 
tions in (47) and (48). In the following, we derive simple design con- 
straints such that when these constraints are satisfied, o>„ will satisfy 
the specifications in (47) and (48). 

So far, we have considered the arbitrary starting sequence described 
at the beginning of this section. Since we may always designate the 
station that is started first as Station 1, we need to consider only the 
following starting sequence in the sequel: At an arbitrary time U < 0, 
the counter at Station 1 is activated. The counter at Station 2 
is activated at t = 0. 

There are two cases to be considered: co, ^ w 2 and to, > co 2 . Our 
analyses of these two cases yield the same design constraint; hence, 
we describe only the case o>i ^ o> 2 . 

Note from the starting sequence that for t s; 0, Station 2 is the 
master and Station 1 is the slave. We therefore can use the results in 
Section III to bound »,(<) for t ^ 0. From this, we can show that oj 
always satisfies the following inequalities: 



COo > 



, (fci«i + fcia 1 fc 2 a 2 T 21 )(co2 — coQ 
1 kycti + k 2 a 2 + k l aik 2 a 2 (T 12 + r 2i ) 

_ [ffc 2 a 2 + k 1 a l k 2 ot 2 T l2 ]§e i ai 

fcxttj + k 2 a 2 + /c 1 a 1 fc 2 a 2 (Ti2 + t 2 ,)J 



and 



too < w 2 + 



[|fc 2 q 2 + fc 1 ai/c2a 2 Ti 2 ]6eia l 
kiot, + k 2 a 2 -\- A; 1 a 1 fc 2 a 2 (ri2 + r 2I ) 



(49) 



(50) 



It should be clear from (49) and (50) that, regardless of the values 
of the time delays r 12 and r 2i , one can easily select the gain e^i of 
the first station so that cu„ will satisfy the constraint in (47). To show 
this more explicitly, we further simplify (49) and (50) (this simplifica- 
tion will, however, make the constraint on e,ai slightly more stringent). 
Since w satisfies (49) and (50), w will definitely lie in the following 
broader range 

o>i — 6e,ai < o> < w 2 + 6eiai . (51) 

Comparing (51) with (47) shows that u satisfies the specification 
in (47) if 



e ' ai < 6ft 



(52) 



1900 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1972 

From (52), one can easily determine the value of e 1 a 1 . Since we have 
designated the station that is started first as Station 1, and since either 
station can be started first, (52) should be applied to both stations. 
To emphasize this, we replace (52) with the following two constraints 

M < £ (53) 

e 2 a 2 < ^- (54) 

Now, to summarize this section: we have shown that if the gains 
of the two stations are designed to satisfy the simple constraints in 
(53) and (54), w will satisfy (47) and the steady-state signaling rates 
will be satisfactory. Since (53) and (54) can be easily satisfied, and 
are independent of all the other parameters in (46), we conclude that 
the steady-state signaling rate can be easily made satisfactory re- 
gardless of the starting sequence, the initial system conditions, and 
the time delays in the communication channels. 

V. SUMMARY AND CONCLUSIONS 

Sections II and III examine the behavior of the system in the master- 
to-slave mode. The station with the slave clock (Station 1 in Fig. 1) 
represents a data station, while the station with the master clock 
(Station 2) represents a station in the digital transmission facility. 
The slave clock at Station 1 employs a digital rate-locked loop which 
consists of a digital counter, an RC filter, and a slave oscillator (Fig. 2). 
The counter is not restricted to have only one stage. A mathematical 
model of the system is formulated in Section II. Transient response 
of the system is determined in Section 3.1. It is shown that, under 
the condition 1 — ICRkcti > in (12), the signaling rate of Station 1 
approaches that of Station 2 in a monotone fashion (transients either 
decay exponentially as shown in (13), (14), (15) and (16), or can be 
tightly bounded as shown in (19)). 

From the transient response, settling time of the slave oscillator 
can be easily estimated. As discussed at the end of Section 3.1, this 
settling time can be rather long when the RC filter has a large time 
constant. For fast start-up purpose, it may be desirable for Station 1 
to transmit data before the slave oscillator is completely settled. Thus, 
during the start-up period, data can be transmitted from Station 1 
to Station 2 faster than it can be transmitted out of Station 2. Con- 
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sequently, a buffer storage is required at Station 2. This buffer size is 
determined and is given in (25). Section 3.4 examines the size of the 
counter in the rate-locked loop (counter size determines the pull-in 
range of the rate-locked loop). In order to avoid counter overflow 
(that is, to ensure pulling in), the counter must have a certain minimum 
size. This minimum size is determined and is given in (26). 

As emphasized in Section 3.2, the slave oscillator in the rate-locked 
loop locks to neither the instantaneous frequency nor the phase of 
the master oscillator. It locks only to the rate of zero crossings of the 
master oscillator. For this reason, we refer to this control loop as a 
rate-locked loop, instead of a frequency-locked loop or a phase-locked 
loop. This difference, while immaterial in the present application, 
should be carefully noted in other applications. 

Section IV examines the behavior of the system in the slave-to-slave 
mode. The two stations to be mutually synchronized represent two 
data stations connected by analog transmission facilities. A rate- 
locked loop is used at each station, and an RC filter is included in 
each loop. A random starting sequence is considered where either 
station can be started first, with the other station activated at an 
arbitrary later time. When the two stations are mutually synchronized, 
the two stations settle to the same steady-state signaling rate hu a 
(h is a proportionality constant and a> is given in (46)). Equation 
(46) shows that o>„ depends on the gains of the counters and oscillators, 
the initial conditions of the counters, filters, and oscillators, and the 
time delays in the communication channels. It is shown that, although 
o) depends on so many parameters, the steady-state signaling rate 
ho3 will lie within desired limits if the simple design constraints in (53) 
and (54) are satisfied (these conditions can be relaxed by using the 
more complicated equations (49) and (50)). These results show that 
the steady-state signaling rate of the system can easily be made satis- 
factory regardless of the starting sequence, the initial system conditions, 
and the time delays in the communication channels. Therefore, there 
is no need to attempt to activate the two stations simultaneously or 
to equalize the delays and gains of the communication channels. 

In conclusion, the detailed transient and steady-state analyses 
show that a synchronization system employing digital rate-locked 
loops can be designed to operate successfully both in the master-to- 
slave mode and in the slave-to-slave mode. Such a synchronization 
system, therefore, is useful in applications where both digital and 
analog transmission facilities are utilized in connecting data stations 
or other types of terminals. 
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APPENDIX 

In this appendix, we first introduce the concept of equilibrium. 
The system is said to be in equilibrium if, corresponding to every 
digit received from Station 2 (the station with master clock), Station 1 
(the station with slave clock) also transmits a digit back to Station 2. 
Then we prove a general theorem which states that the system will 
reach equilibrium if the arbitrary filter F(s) (not necessarily an RC 
filter) satisfies a simple condition. Based on this general theorem, we 
then show that when an RC filter is used, the signaling interval of 
Station 1 will lock to that of Station 2 exactly. 

For brevity, we define p,(0 as co t t + a, j Vi{r)d,T -f- 0i . The lowpass 
filter transfer function F(s) can always be normalized such that F(0) = 1. 
Clearly, any useful lowpass filter must cut off as frequency approaches 
infinity; therefore, we can write F(<x>) = 0. By changing units, we 
can and shall set e x = ir and ai = 1. Without loss of generality, we 
assume that co 2 — cj, > 0, and that the counter counts both upward 
and downward zero crossings. The zero crossings of Sj(£) and s 2 (t) 
control the signaling rate of Station 1 and Station 2, respectively. 
Let T be the time interval between each two consecutive zero crossings 
of s 2 (t), that is, T = tt/w 2 . When the time interval between each two 
consecutive zero crossings of Si(t) also becomes T, signaling rate of 
Station 1 locks to that of Station 2. Thus, to determine the locking 
behavior, we need only to examine N t (t) when t — > oo. Since N t (t) 
can be deduced from p t (t), v t (t), or Ui(t), we shall examine either pi(t), 
or v t (t), or Ui(t) in the following analysis (depending on which one is 
the most convenient). 

The behavior of the system is governed by the equation 

P.(0 = Wi* + f [/*w,(p,)] clr + 0, (55) 

where * denotes convolution, and the symbol Mi(pi) indicates that Ui 
is a function of p, . Since w x depends on p, through the nonlinear zero- 
crossing counting process, (55) is a nonlinear differential-integral 
equation. It is impossible to solve this equation for all t, so we shall 
first examine v,(t) and Ui(t) to obtain a steady-state solution of this 
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equation. Then we shall consider the uniqueness of this steady-state 
solution. 

From the mathematical formulation in text, 



(56) 



Mx(«) - (o> 2 - Wi)< - / Vt(r) dr + B 2 - 6, + *,(*) - * a (oJ ■ 

Let Uiis), V 1 (s), ¥,(s), and ^ 2 (s) be the Laplace transforms of m x (0, 
»i(0i ^i(0, and ^,(0, respectively. From (56) and 7,(s) = F(s)U 1 (s), 
we obtain 

*« - #nfe + £n=s + r6» *' (s) " ^> ** (s) (57) 

and 

v M F(s)(co 2 - co.) F(s)(fl 2 - 0.) , F(s> ^ , , _ g@g * () 

(58) 

We wish to determine the Ni(t) that satisfies the system equation 
(55) when t — > °o . For brevity, such a solution is called a steady-state 
solution. From (58), a steady-state solution is obtained, and is stated 
in the following theorem. 

Theorem 1: At steady-state (that is, when t —> oo), (55) is satisfied if 

NS) = #■(« - r„) (59) 

where t is such that the mean value oj u x (t) is w 2 — «i • 

Proof: Since fa(t) and ^ 2 (0 do not approach a limit when t — ► », 
one cannot apply final value theorem to the last two terms in (58). 
However, final value theorem can be applied to the first two terms. 
This yields 

*« - - - - + H>?fe*'«] - **\j+h™] ■ 

t-* oo. (60) 

The condition / -» oo applies to the rest of the proof. Clearly, (59) 
is equivalent to the statement that 

P,(0 = <*A + p,(0 (61) 

where p\(0 is a periodic function of period T. Thus, to prove Theorem 1, 
one needs only to show that the right side of (55), when computed 
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from (59), is identical with the right side of (61). From (59), fa(t) is a 
periodic function of period T. Consequently, £~ l {[F(s)s/s + F(8)]¥i(a)} 
is a zero-mean periodic function of period T. Since fa(t) is a period 
function with period T, £~ 1 {[F(s)s/s + F(s)\9 a (s)} is also a zero- 
mean periodic function of period T. Thus, from (60) 

Vi(t) = oj 2 - wi + »i(<) (62) 

where Si (t) is a zero-mean periodic function with period T. Substituting 
the Vi(t) in (62) for the integrand [/ * Wi(pi)] in (55), we see that the 
right side of (55) is identical with the right side of (61). This proves 
Theorem 1. 

Equation (59) in Theorem 1 implies that signaling rate of Station 1 
locks to that of Station 2. Now we consider the problem of uniqueness 
(that is, whether (59) is the only steady-state solution). We first prove 
that, under a simple condition, Station 1 cannot add or delete bits 
from a customer's data stream. 

As described in Section II, the zero crossings of s t (t) and s 2 (t) control 
the signaling rates of Station 1 and Station 2, respectively. More 
specifically, Station 2 transmits the with digit to Station 1 at the raw th 
zero crossing of s 2 (t); and Station 1 transmits the rath digit to Station 2 
at the ran th zero crossing of s t (t) (in practice, n ^ 1). Thus, Station 2 
transmits a digit to Station 1 every n T seconds. We say that the 
system is in equilibrium if, corresponding to every digit received from 
Station 2, Station 1 also transmits a digit back to Station 2. More 
precisely, the system is in equilibrium if we can partition the time 
axis into A^-second time intervals such that Station 1 will transmit 
a digit back to Station 2 in each of the n T-second time intervals. 
Theorem 2: The system ivill reach equilibrium ij 

-' < *i+m* M ] K '■ (63) 

Proof: The condition / — > oo is implied throughout this proof. Using 
the final-value theorem, one can show from (57) that when t — > oo, 

«i0) - « 2 - «1 - T 2 (t) + (7,(0 (64) 

where 

' M = £ "'k+Tw*' (s) ] (65) 

°' (l) - £ ~'k+iw *■<«>] (66) 

and £ _1 denotes inverse Laplace transform. 
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Since \{/ 2 (t) is periodic with period T, <x 2 (t) is a zero-mean periodic 
function of period T. Let max cr 2 (0 and min o 2 (t) be the maximum 
and minimum value of o 2 (t), respectively. We now determine max o 2 (t) — 
min o- 2 (0- Note that yf/ 2 (t) can be written as 

UO = M + e 2 ) - £ «*[« - («' + «01 ( 67 ) 

where u(t) is the unit step function denned by 

m(0 = 0, < < 

= 1, t > 0. (68) 

When ^ 2 (0 is applied to a network with transfer function [s/s + ^(s)] 
(hereafter called network A), the output is a 2 {t). Clearly, when the 
first term u 2 t + 62 in (67) is applied to network A, the output is a 
continuous time function for t > 0. The second term in (67) consists 
of unit step functions. It can be shown that, when a unit step function 
u(t) is applied to network A, the output is unity when t = + , ap- 
proaches zero when £—><», and is continuous f or < / < °° . From 
these results, it is clear that 

max o- 2 (0 — min a 2 (t) ^ it. (69) 

We have set e x = t. Therefore, Ui(t) is a multiple of ir. We are con- 
sidering the case w 2 — w, > 0. As illustrated in Fig. 7, let n be an integer 
such that 

hit ^ o> 2 — u>i < nir + ir. (70) 

It is clear from (70) and (69) that there is a t at which oj 2 — cu, — a 2 {t) 
equals nir or (n + l)ir (let this / be denoted by t t ). Note that o> 2 — oil — 
o- 2 (0 rnay intersect only the level nir, or only the level (n + l)7r, or 
both the levels. For this proof, we need to consider only the first case. 
Since a 2 (t) is periodic with period T, u 2 — an — o- 2 (0 is also periodic 



(n+ 1)7T 



(D 2 - O} 



t!-T t, ti+T ti + 2T 

Fig. 7 — Illustration for the proof of Theorem 2 (showing the definition of n and 
the partition of the time axis into successive T-second intervals). 
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with period T. Thus, co 2 — wi — tr 2 (t) must intersect the level nir also 
at time instants U + iT, i = ±1, ±2, • • • . These intersections are 
illustrated in Fig. 7. 

Now consider the value of u x {£) at an intersection t r + zT, i = 0, 
±1, ±2, • • • . From (63), we have, at t = t t + »T 

fl* — 7T < Mi(0 < ??7T + 7T. (71) 

Since w,(0 must be a multiple of 71- , (71) implies that 

U&) = m? (72) 

at t = <i + iT, t = 0, ±1, ±2, • • • . Since iV 2 (0 increases by one every 
T seconds, (72) requires that N x (t) increase by one in each of the T- 
second intervals illustrated in Fig. 7. This proves Theorem 2. 

Now we consider the case where the filter F(s) is the usual RC filter. 
We first prove that eq. (63) is satisfied in this case (consequently, the 
system will reach equilibrium). 

When RC filter is used, 



where 



7+Txij »•« = *'< s > - CR( S + ,!x. + rj *' (S) (74) 



- - 1 + A £7^ W» 



1 - Vl - 4Cfi 
r2 = 2CR (76) 

The system is designed such that 

1 - 4CR > 0. (77) 

Therefore ? - i and r 2 are real numbers and 

r, > r 2 > 0. (78) 

Let [l/CR(s + ri)(s + r 2 )] be denoted by G'(s), then 

g(t) = £-\G(8)] = c/g(ri 1 _ r>) [e- r " - <T r "]. (79) 

From (78) and (79), 

C(0 > 0, t > 0. (80) 
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(81) 



(82) 



(83) 



From (80), we can write 

= [ tMq(t - t) (It < I rrgit - t) (It < / Tf/(r) dr. 

Jo J o J ° 

Clearly, J" (j(t)(It = 6'(0) = 1. From this and (81), we have 

< £~*[*i(«)G(«)] < *■ 
From (82) and (74) 

Hence, (63) is satisfied and the system will reach equilibrium. 

Next, we examine the detailed behavior of the rate-locked loop. 
Note that there are two basic variables in the rate-locked loop, namely, 
Ui(t) and w,(0- Let the Wi(0 and Vi(t) corresponding to the steady-state 
solution in (59) be denoted by «?(*) and v*(t), respectively. First, we 
sketch u*i(0 and v\(t). From Section II in text, s 2 (t) = sin [<a a t + 2 ]. 
To simplify our graphs, let us omit 2 . Then N 2 (t) jumps by 1 at t = 
IT, I = 0, 1, 2, • • • . From this and (59), we see that w?(0 is as sketched 
in Fig. 8, where I denotes an arbitrary integer. The pulse width y* 
in Fig. 8 is such that the mean-value of u1(t) is w a — o>i . Therefore, 



T 

y* = — [co 2 — co, — nir]. 

TV 



(S4) 



u,*(t) 
/ 



(n+ 1)7T- 



O) 2 -a)i 



- u" 



t*-v-^i 



— 1>* 

!&>. — 



K-y-»l 

I I 



-■"V 



-v,'(t) 



£T 4T+T 

Fig. 8— Sketch of «i*(0 and v t *(l). 
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Since u*i(t) is periodic with period T, v\(t) is also periodic with period T. 
As can be seen from v\{t) in Fig. 8, u*(t) charges the capacitor C in 
the time interval IT to IT + y*, and the capacitor C discharges in the 
time interval IT + y* to IT + T. Let v* denote the value of vf(0 
at t = IT + y*. Clearly, v* must have such a value that v*(t) has a 
mean-value of ta% — <a x . 

In order to show that u\{t) and v\(t) are the only steady-state solu- 
tion, we begin by assuming different Ui(t) and v y {t), and demonstrate 
that they must approach u\(t) and v\(t) as t increases. We have proved 
that the system must reach equilibrium. From (72), when the system 
reaches equilibrium, Ui(t) = nv at t — ti + iT, i = 0, ±1, ±2, • • • . 
[As can be seen from the discussion after (70), Ui(t) may assume the 
other value (n + l)ir at such time instants. However, these two cases 
are similar and we need to consider only the first case.] Therefore, Ui(t) 
can assume only one of the two forms in Fig. 9 in each of the time 
intervals t t + iT to U + iT + T. The first form is illustrated in the 
time interval t\ to t% + T in Fig. 9, while the second form is illustrated 
in the time interval t x + T to U + 2T. In the first form, the zero crossing 
of Si(t) (represented by the downward arrow) takes place prior to the 
zero crossing of s 2 (t) (represented by the upward arrow). The order 
is reversed in the second form. Note that, if m x (<) always assumes the 
first form, one would have v t (t) < nir. From this, one can easily show 
that Ui(t) cannot always assume the first form in the successive T- 
second intervals. Next, consider the width of the pulse when Ui(t) 



(n+1)7T 



(n-1)7T 









*-v— * 
















# 




1 1 
1 1 
1 1 
1 1 
1 1 


1 1 

1 1 1 

1 1 1 



\\ ti+T tj+2T 

Fig. 9 — Illustration of the two forms of Ui(t). 



DATA-RATE SYNCHRONIZATION SYSTEM 



1909 
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Fig. 10— Sketch of ui(t) and vi(t). 



assumes the second form. This width, designated by y in Fig. 9, may 
vary from one T-second interval to the next. If this width were always 
less than y*, v^t) would be less than v* for all t. From this, one can show 
that this width cannot always be less than y*. From these results, 
there must be some T-second intervals in which Ui(t) assumes the 
second form and the pulse width y is equal to or greater than y*. We 
shall select one such time interval (say, the time interval h + iT ^ 
t < t { + iT + T illustrated in Fig. 10) and examine u x (t) and v^t) 
for t > ti + iT. We need to consider only two cases (refer to Fig. 10): 

Case 1: v l (t l ) < v* 

Case 2: v l (t l ) ^ v*. 

The instantaneous radian frequency of VCOi is a>i + aiV^t), where 
w, is the free-running radian frequency and a^t) is the correction 
term. In data communications, o^ is very close to the radian frequency 
o) 2 of the master clock. (For example, it may be specified that the 
maximum difference between w, and w 2 be limited to 0.005 percent 
of fa> 2 .) Consequently, only a very small correction term a r v v {t) is needed. 
For this season, the time interval between each two consecutive zero 
crossings of s t (t) is essentially determined by the term a>i£ in p t (t). 
Therefore, the pulse width y changes only very slightly from one pulse 
to the next (in other words, in Fig. I0y t+ i is very close to y t ). 

For the purpose of illustration, in Fig. 10 v^t) is shown to increase 
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and decrease quite rapidly in each T-second interval. In practice, the 
filter time constant RC is several orders larger than the time interval T 
(for example, RC = 10" ' seconds, T = 10~ 5 seconds). Thus, v t (t) is 
essentially a constant in each T-second time interval. 

Now consider y l + i and v^ti+i), 3 = 1> 2, 3, • • • .It can be shown 
rigorously that if there is an h such that 

yi + h - y* 

*>i(Jj + a) = v* 

then y t+i = y* and Vi(ti + i ) = v* for all j > h. Therefore, to show that 
Uiit) and Vi(t) approach W\(t) and V\(t), we need only show that y l+i 
and Vi(ti +i ) approach y* and v*, respectively. 

Now consider Case 1; after ti + %T, yi + j and Vi(t t + i ) approach y* 
and v* in three stages. Immediately after ti + iT, Vi(t l+ j) is less than 
v*. Consequently, the time interval between each two consecutive 
zero crossings of Si(t) is slightly larger than T, and y i+i increases slowly 
with j. (Note from Theorem 2 that y l + i must remain less than T.) 
Since Wi(<i + ,-) is less than v* and y l+j remains larger than y*, Vi(t t + j) 
must increase slowly with ;. The second stage starts when v x (ti + i ) 
reaches u*. Since the pulse width y l + i is larger than y*, Vx{ti +i ) keeps 
increasing with j. (Note from Theorem 2 that Vj(t l + i ) cannot exceed 
(n + l)r.) This implies that Vi(ti + i ) will be larger than v*. Consequently, 
y l+i must decrease with j. Clearly, when y i+i decreases, the rate of 
increase of t>i($i+,) decreases. The third stage starts when yi + i de- 
creases to such a value (still larger than y*) that v 1 (i J+J ) ceases in- 
creasing. Since Vi(ti +i ) is now larger than v*, y l + i must keep decreasing. 
Clearly, this must also reduce t>i(tfj +J -). Consequently, y l+i and Vi(t t+i ) 
approach y* and v*, respectively. 

The above discussion is for Case 1. It can easily be extended to 
Case 2. Thus, when the system reaches equilibrium, the time interval 
between each two consecutive zero crossings of Si(t) will be exactly T 
seconds. 
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